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Abstract 

A general spin symmetry argument is proposed for spin currents in semiconductors. In particular, due to 
the symmetry with respect to spin polarization of the helicity eigenstates of the Luttinger Hamiltonian for a 
hole-doped semiconductor, the spin polarized flux from a single helicity eigenstate induced by an external 
electric field, is canceled exactly when all such contributions from eigenstates that are degenerate in energy 
are summed. Thus, the net spin current predicted by Murakami et al, Science 301, 1348 (2003), cannot 
be produced by such a Hamiltonian. Possible symmetry breaking mechanisms which may generate a spin 
current are discussed. 
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Although generating and utilizing spin currents in metallic films and tunnel junctions has be- 
come a routine matter in spintronics research and application, generating spin currents in semi- 
conductors has attracted greater research efforts because of the tremendous potential of combining 
spintronics with the well-developed semiconductor technology. However, lack of an effective 
method for spin injection into semiconductors has been the main obstacle of progess. The recent 
work by Murakami et al 1 is remarkable that it points an unexpected possible route for achieving 
spin injection. They proposed that a spontaneous and dissipationless spin current would be pro- 
duced by an applied electric field in a hole-doped semiconductor, thus bypassing the difficulty 
of traditional methods of spin injection. Other groups 2 are considering alternative proofs for the 
existence of such a spin current. The unconventional and extraordinary nature of this idea, and 
its potentially wide impact in the field of condensed matter physics, calls for a careful and skep- 
tical scrutiny. In this letter, we will examine the possibility of a spontaneous spin current from 
the consideration of the spin symmetry of the electron wave function. Symmetry rules are more 
general than the particular form of the Hamiltonian one solves. Although our discussion will be 
focused on the Luttinger Hamiltonian, for which a negative result is obtained, the approach is 
equally applicable to other forms of Hamiltonian that might be considered for spin transport. 

The result of Murakami et al was based on the consideration of the Luttinger Hamiltonian 3 , 



which describes the valence bands for a large class of hole-doped semiconductors. Note here that 
S represents the total angular momentum of an electron. The expectation value of the spin part of 
the angular momentum is related to S through (l/3)(ip\S\ip). Under the influence of a uniform 
electric field, they obtained a result that a net spin current is induced in a direction perpendicular 
to the field. This is a puzzling result because the size of the spin current seems to be independent 
of the strength of the coupling that is the origin of this spin current, the coefficient 72 in Eq. ©. In 
fact, the spin current becomes a constant when 72 approaches zero which seems to predict a spin 
current in an isotropic free-electron gas. 

The key to the solution in Ref. 1 is the 4x4 unitary transformation U (k) which generates the 
helicity eigenstates, 



by the diagonalization U (k) (k • S)f7(k)t = kS z . Here Si, i — x,y, z, are the 4 x 4 spinor matrices 
for spin 3/2 (not to be confused with the zth component of the spin matrix S). Although the choice 
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(1) 



(k- S)|k, A) = fcA|k,A), 



(2) 



of U(k) is not unique, Ref. 1 chose the simple form U(k) = exp (i9S y ) exp (i(f)S z ). This represents 
a rotation in the spin space of angle <f> about the z axis followed by a rotation of angle 6 about the 
y axis. Under this transformation, the coordinate system for the spin part of the wave function 
has been rotated to a reference frame which is local in the k space, and the good quantum number 
is the helicity A within this frame. The A eigenstate corresponds to a spin state that is precessing 
around an axis parallel to k, and A is the projection of the spin along this axis. The ±A eigenstates 
at the same k point have exactly opposite spin polarizations. The spin symmetry of the helicity 
eigenstates are expressed as, 

(k, A|5<|k, A) = -<k, -A|5i|k, -A), (3) 

where % = x, y, z. A time-dependent spin operator, S(t) = S x (t)x. + S y (t)y + S z (t)z, can always 
be expressed in the form with the ith component, Si(t) = J2j ®ij(t)<Sj, where oty (t) are elements 
of a rotation matrix and are scalar functions of time. Thus, the antisymmetry in A holds for any 
time-dependent spin operator, 

(k, A|S(£)|k, A) = -(k, -A|S(£)|k, -A). (4) 

Because the ±A states are degenerate in energy, the net spin from each k direction is exactly zero. 
This is consistent with the fact that semiconductors not doped with magnetic impurities generally 
do not exhibit magnetism. Similarly, using the following symmetrized three-spin product which is 
also antisymmetric in A, 

(k, \\{S h {S„ S k }}\k, A) = -(k, -\\{Si, {S„ S k }}\k, -A), (5) 

with i, j, k = x,y, z, we find the symmetry relationship, 

(k, \\{Si(jt), {Sj(jt), S(t)}}|k, A) = -(k, -\\{Si(jt), {S 3 (t), S(t)}}|k, -A), (6) 

which will be needed in computing the spin current. 

The question remains whether the spin polarization is canceled in a linear response solution 
under a uniform electric field. Let us consider the transverse component (perpendicular to the 
applied field) of the linear-response current which arises from the first term on the rhs of the 
equation of motion, 
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where Xi = U(k)xiW(\c) and we use a tilde to clearly distinguish between the transformed op- 
erators Xi and the untransformed operators Xi. We will consider the meaning of the second term 
later in this paper. Karplus and Luttinger 4 pointed out in their consideration of the anomalous Hall 
effect, that a transverse current is the consequence of the spin-orbit coupling which introduces a 
left-right asymmetry in the stationary (Bloch) states of the system. The transverse current here 
similarly arises from the (k • S) 2 term in the Hamiltonian. What is remarkable about the transverse 
current predicted by Ref. 1 is that for each direction k the contribution to the transverse current is 
completely independent of the parameter 72, the coefficient of the (k • S) 2 term in the Hamilto- 
nian. Thus it makes an unusual prediction of a finite spin current in the direction perpendicular to 
the external electric field even for an infinitesimal coupling constant 72 . If this were indeed true, 
it would cast doubt on the validity of the linear-response theory itself. Fortunately, because the 
equation of motion is the same for both ±A eigenstates, the transverse currents for both states are 
the same. Since the spin polarizations of the two states are exactly the opposite, there is no net 
spin current arising from this transverse term. 

The total spin current can be calculated directly within the unrotated global reference frame. 
The equations of motion without the rotation operations are simply 

k(t) = eE, (8) 

and 

±(t) = — (71 + -72) - -7 2 [(k • S)S + S(k • S)]. (9) 

m 2 m 

The electron flux at each k is Z^(k, A|x|k, A). When summed over all k, the total flux is cancelled 
except for the drift contribution due to the combination of Eq. ® and scattering, which within the 
linear response theory yields a term proportional to the electric field E. If we set the z direction as 
the direction of the applied electric field, then the transverse spin flux can be calculated using the 
helicity eigenstates, |k, A), of (k ■ S) which are also the eigenstates of the Hamiltonian in Eq. ([T]). 
The spin current is calculated from the linear response formula Tr[pj spin ] where p is the density 
matrix and the spin current operator is given in the symmetrized form, 

jspin = ^{S,X}. (10) 

This would be the typical starting point for deriving the Kubo-Greenwood linear response for- 
mula. However, in our case it is more straightforward to use the Heisenberg representation, and 



4 



directly apply the equations of motion, Eqs. dU) and ©. The flux in the y direction with the spin 
polarization along x at each k point, is 

fyiV = ^E<M|{^(^(t)}|M>n A (k) 

A 

= |^(7i + 572) £<k,A|^)|M)n A (k) 

A 

-^7 2 k • £(k, A|{^(t), {^(t), S(t)}}|k, A)n A (k), (11) 
6m ^ 

where S^t) and S y (t) are the components of the spin matrix in the x and y directions, and n x (k) is 
the occupation number for the state |k, A). Before summing over A, a linear response spin current 
would arise from the drift term in k. However, the spin factors can be shown to cancel exactly at 
each k. Because of the degeneracy of the helicity eigenstates |k, A) and |k, — A), the occupation 
number, n x (k) = n~ x (k) is symmetric with respect to A. Using the symmetry relations, Eqs. © 
and ©, we obtain, 

= 0. (12) 

Thus no transverse spin flux is induced by applying an electric field. 

Now let us consider the second term in Eq. Q. We show here that it results from the time 
derivative of the local reference frame. First, let us rotate the equation of motion for x, Eq. ©, 
to the local reference frame. In order to arrive at an expression comparable to Eq. ©, one needs 
to invoke the adiabatic approximation to drop the off-diagonal terms coupling the light-hole (LH) 
and heavy-hole (HH) bands in U(k)SiW(k) = hS-Jk. We find, 

C/(k)^C/t(k) = ^( 7l + | 72 - 2Sh 2 ) = ^. (13) 
m 2 nix 

This is exactly the first term in Eq. ©. It is clear, then, that the time derivative of x, includes two 
contributions, 

it = U(k) Xi U\k) + k • [V k f/(k)a; l f/ t (k) + U(k) Xl V^ (k)], (14) 

where the second term is equivalent to the second term in Eq. © and results from taking the time 
derivative of the rotation matrix U . This clearly shows that x and x are not related by a gauge 
transformation, the latter including a contribution from the time derivative of the reference frame. 
Therefore it would be incorrect to compute the spin current from x as is done in Ref 1 . In fact, the 
validity of the adiabatic approximation is also questionable. However, this approximation is only 
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needed to clearly separate x into two terms. It remains true that x and x are not related by a gauge 
transformation whether or not the adiabatic approximation is invoked. We thus conclude that the 
nonzero spin current is an artifact of their approach, not a physical solution. 

A further comment on the generality of our result is in order. We have shown that the cancel- 
lation of the spin current is due to the degeneracy of the ±A helicity eigenstates of the Luttinger 
Hamiltonian, Eq. CO). One is certainly tempted to consider the possibility that by breaking the 
degeneracy between the ±A states, one might be able to produce a spin current from Eq. (fTTT) . For 
example, one might introduce an anisotropic kinetic energy term in Eq. ([T]), such that it breaks the 
spatial inversion symmetry, thus Eq. (fTTl) will need to include interband terms which are usually 
not symmetric with respect to ±A. Another example that breaks the symmetry between the ±A 
states has been considered by Ref. 2 , using the Rashba Hamiltonian which is linear in k and S. The 
negative result for the Luttinger Hamiltonian here does not rule out the possibility of a spin current 
when such symmetry is broken. 

In conclusion, we have shown that at least for a class of semiconductors described by the 
Luttinger Hamiltonian, spin symmetry of the eigenstates rules out the possibility of a spontaneous 
spin current in these materials. Thus, any attempt to produce such a spin current must include 
symmetry breaking terms in the Hamiltonian. This should provide important guidance to future 
attempts of finding a mechanism for spin currents in semiconductors. 
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